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ABSTRACT 

We predict the flux and surface velocity perturbations produced by convectively 
excited gravity modes (g-modes) in main sequence stars. Core convection in massive 
stars can excite g-modes to sufficient amplitudes to be detectable with high precision 
photometry by Kepler and CoRoT, if the thickness of the convective overshoot region 
is < 30 per cent of a pressure scale height. The g-modes manifest as excess photo- 
metric variability, with amplitudes of ~ 10 micromagnitudes at frequencies 10 /zHz 
(0.8 d~^) near the solar metallicity zero-age main sequence. The flux variations are 
largest for stars with M > 5 M©, but are potentially detectable down to M 2 — 3 
Mq . During the main sequence evolution, radiative damping decreases such that ever 
lower frequency modes reach the stellar surface and flux perturbations reach up to ~ 
100 micromagnitudes at the terminal-age main sequence. Using the same convective 
excitation model, we confirm previous predictions that solar g-modes produce surface 
velocity perturbations of < 0.3mms~^. This implies that stochastically excited g- 
modes are more easily detectable in the photometry of massive main sequence stars 
than in the Sun. 
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1 INTRODUCTION 

At the interfaces between convective and radiative zones 
in stellar interiors, convective motions transfer a fraction 
of th eir kinetic energ y into wa ves in the radiati ve layer 
(e.g., [Pres s 1981; Gol dreich fc Ku mar 1990; Bclkac em et al.l 
|2008| ). Studying the details of this energy transfer is impor- 
tant for our understanding of mixing at convective bound- 
aries, the evolution of shear layers and the excitation of stel- 
lar oscillations in many contexts. 

The Sun provides an exquisite laboratory for studying 
convectively excited sound waves (p- mod es), with a forest of 
modes observed at the solar surface (see lGizon et al.|[2010l : 
IChristensen-DalsgaardI |2002| . for recent reviews). These 
sound waves have been detected in many other stars as well, 
with rapid growth in the sample of observed "solar-like os- 
cillators" in the past several years thanks to high precision 
monitoring campaigns like Kepler (|Beddinell201ll ). 

To date, gravity modes (g-modes) have not 
been convincingly obser ved at the solar surface (see 
lAppourchaux et al.l I2OI0I . for a recent review). Gravity 
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modes are observed in several other stellar types, inclu ding, 
for example, white dwar fs (jWinget fc Kepleil |2008| ) and 
slowly-pulsating B-stars (|De Cat! 20071 ). In each of these 
cases, the g-modes are linearly unstable, rather than 
stochastically excited as in the case of solar- like oscillations. 
In the interesting case of the mixed gravit y and pressure 
modes recently observ ed in giants (see e.g., iBedding et al.l 
l201ll : lBeck" et a l. 2011), the convective excitation occurs in 
the envelope where the modes behave locally as pressure 
waves. Thus, the direct excitation of g-modes by turbulent 
convection is not well tested observationally. 

Main sequence stars more massive than the Sun have 
convective cores and relatively compact, predominantly 
radiative envelopes, through which g-modes propagate. 
This structure makes them potential hosts for observable, 
stochastically excited g-modes. The unprecedented micro- 
magnitude precision of photometric monitoring campaigns 
like Kepler (Koch et al.i.2010, ) a nd CoRoT (Convection , Ro- 
tation and planetary Transits; lAuvergne et al.l |2009| ) pro- 
vides exciting prospects for detecting these stochastically 
excited g-modes and opening a new window into both the 
physics of convect ive bound aries and massive stellar inte- 
riors. Indeed, iSamadi et al.l (|2010i ) studied the convective 
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excitation of g-modes in main sequence stars with masses of 
10, 15, and 20 M© and central hydrogen mass fraction of 0.5, 
finding that convectively excited g-modes in these stars may 
reach amplitudes near the threshold for detectabihty with 
CoRoT. In the following, we present a complementary ap- 
proach and investigate g-mode amplitudes in a wider range 
of initial stellar masses, from 2 to 30 Mq , and evolutionary 
states along the main sequence. 

We begin by describing our model for the spectrum of 
g-mode excitation in 321 followed by our method for calcu- 
lating the observable signatures of convectively excited g- 
modes in ij3] [l] and [5] In ijG] we present the results of our 
calculations for the surface flux and velocity perturbations of 
g-modes in massive main sequence stars. We conclude with 
a discussion of remaining uncertainties, including the effects 
of rotation and stellar evolution modeling, and the current 
and ongoing observations which might shed light on these 
predictions (®. 



2 CONVECTIVE EXCITATION OF GRAVITY 
MODES 

Convection efficiently excites gravity modes in an ad- 
jacent stably stratified medium, as both convective 
and g-mode pert urbations are roughly incompressible 
l|Goldreich fc Kumar 1990 1. The power spectrum of g-modes 
excited by stellar convection at and above the characteristic 
convective turnover frequency, ujc, has been calculated us- 
ing both heuristic physical arguments (setting the wave pres- 
sure in the stable layer e qual to the convective ram pressure; 
|Presslll98ll : iGarcia Lopez fc Spruit. 19911 ) and by solving the 
inhomogene o us wave equation with convective source terms 
|Presdll98ll; 'Goldreic h fc Kumaij [l990l : iKumar et all 1 19991 : 
[jBelkaccm et al. 2009a). These approaches are reviewed and 
ex tended to more realis t ic rad iative-convective boundaries 
in iLecoanet fc QuataertI (|2012l ) ; we summarize here the key 
results for our application. 

The low frequency g-modes that carry most of the 
wave power are strongly damped by radiative diffusion (see 

As a result only g-modes with frequencies significantly 
larger than the characteristic convective turnover frequency 
(i.e., uj S> ujc) can set up global standing waves poten- 
tially detectable at the stellar surface - this is true for both 
solar-type stars with envelope convection zones and massive 
stars with core convection zones. High frequency g-modes 
have characteristic radial wavelengths in the radiative zone 
that are large compared to the thickness of the radiative- 
convective boundary. In this limit, it is reasonable to approx- 
imate the radiative-convective bounda ry as a discontinuity 
in th e Brunt- Vaisala frequency (as in I Goldreich fc Kumaij 
Il990t ) . For g-modes excited by a core convection zone (or its 
overshoot region), which is relatively unstratified, the exci- 
tation occurs at a roughly fixed position at the edge of the 
convection zone. The power spectrum of energy supplied to 
g-modes is then given by 
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where Lconv is the total luminosity carried by convection, 
Aic (= I'conv/csound) is the convectlvc Mach number, Vconv is 
the convective velocity, ^max = (r/H) {uj /locY^^ , and r and 
H are evaluated at the edge of the convective core (where 
their ratio is ~ 1). The power-law exponents a and b are 
given by a = 13/2 and 6 = 20 

Note that integrating over frequency and I, equation [1] 
implies that convection exc ites a total g-mode luminosity of 
(|Goldreich fc Kumaj[l990l ) 



Eg ~ Mc 'L. 



c. J-'conv ■ 



(3) 



In fact, the efficiency of g-mode excitation can exceed equa- 
tion [3] depending on the nature of the convective-r adiative 
transition (Press 1 98ll : ILecoanet fc QuataertI [2OI2I ). Equa- 
tion [T] is, however, appropriate for the high frequency g- 
modes of interest in this paper. In our numerical work, we 
choose the uncertain dimensionless constant in equation [T] 
so that the integrated wave power is exactly Eg = Mc Lconv. 

One of the primary uncertainties in the prediction of g- 
mode excitation is whether the excitation occurs primarily 
in the bulk of the convection zone or in a thin "overshoot" 
region near the radiative-convective boundary. We param- 
eterize this uncertainty by taking the local scale-height in 
the excitation region to be rjH, where H is the true pressure 
scale height near the radiative-convective boundary. This 
amounts to taking H — j- rjH in equation [1] when calculating 
the wave excitation. The parameter 77 can take on values in 
the range (0, 1], such that = 1 corresponds to excitation 
dominated by eddies of size ~ If in the convection zone and 
r; < 1 to excitation dominated by a thin overshoot region 
hosting eddies of size ~ r/H. Excitation dominated by a thin 
overshoot region will have a higher characteristic frequency 
iJc/ri since the length scale of the convective motions is re- 
duced by 77. This serves to shift the g-mode power input to 
higher frequency and shorter spatial scales (higher lu and £). 
Note, in particular, that for low £ modes having uu S> i^c, 
Eg oc ff"^''^. Thus excitation in a thin overshoot layer sig- 
nificantly increases the power supplied to the high frequency 
modes that are the most observable. 

In their numerica l simulations of the s o lar ra diative- 
convective transition, ^Rog ers fc Glatzmaieil (|2006l ) found 
that the power into g-modes is roughly constant for uu ~ 
1 — 10 Uc and then decreases significantly for cj > 10 ojc (see 
their figs. 1 fc 2). This increase in the characteristic fre- 
quency of the excited g-modes is reasonably consistent with 
excitation in the overshoot layer given the width of the over- 
shoot layer of ~ .05 H found in the same simulations (see 
iRogers et al.|[200^ . The smaller scales in the overshoot re- 



^ Because excitation of g-modes by Reynold's stresses dominat e 
excitation by entropy fluctuations l lGoIdreich &: Kumar! ll99Cl) , 
the convective luminosity in eqn. [T] is technically that associated 
with the kinetic energy, rather than enthalpy. In mixing length 
theory, these two contributions to the convective luminosity are 
comparable. For the purposes of our calculation, we thus take 
Lconv to be the total convective luminosity in the stellar model 
and absorb uncertainties in the relative contribution of kinetic 
and enthalpy fluxes into the uncertain normalization of eqn[T] 
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gion wou ld also naturally produce pow er at higher I, as also 
found by Ro gers fc GlatzmaieU (|2OO6!)0 

Fig-Eshows the logarithmic excitation spectrum for our 
2 Mq model, for rj — 1 (left panel) and 0.1 (right panel). 
This shows that most of the wave energy is in the lowest 
frequency modes convection can excite. In addition, since 
the wavelength and frequency of convective eddies are cor- 
related, energy input at higher mode frequencies is predom- 
inantly in modes of higher £. 



3 STELLAR MODELS 

We have constructed main sequence models from zero-age 
main sequence (ZAMS) to the terminal-age main sequence 
(TAMS), at a central hydrogen mass fraction of 0.01, for a 
range of initial masses from 2 to 30 using the MES A stel- 
lar evolution code (version 4298; IPaxton et"aLll201lM All 
models are non-rotating a nd solar metallicity (which we take 
to he Z = 0.02) with the iGrevesse fc Sauvali lll998f) chem- 
ical mixture. Following the results of Isrott et aL (|201lh . 
we determine convective boundaries using the Ledoux cri- 
terion with 30 per cent of a pressure scale height of over- 
shoot Q and a mixing length parameter a = 1.5. We also 
use semiconvection with a dimensionless efficiency parame- 
ter, Osc = 0.1, though this is largely irrelevant due to the 
o vershoot. Finally, we assume the th eor etical mass loss rates 
of lVink. de Koter. fc Lamerd (|200ll ) (or lde Jager et al.lll988l 
when Teff < 10* K) scaled down by a factor of 0.8. 

We tested for numerical convergence by varying the 
number of mesh points and timesteps taken during the evo- 
lution, finding good agreement at the level of refinement 
used for the models presented0 We have also tested a range 
of other model parameters, including the mixing length, 
boundary definitions and rotation, to test their effects on 
our results (see ^7.1.3|l . Fig. [2] shows the evolutionary tracks 
in the H-R diagram for the models described above from 
ZAMS to TAMS. Along the sequence, crosses appear every 
1/lOth of the main sequence lifetime, and symbols mark the 
locations where we have calculated the convective excitation 
of g-modes: ZAMS {Xc = 0.68), midMS {Xc = 0.33), and 
TAMS {Xc = 0.01). 

Table Ogives parameters for the convective cores of our 
ZAMS models, including the outer scale convective turnover 
frequency (i^c = cJc/(27r)), the core luminosity (Lconv, which 
is roughly equivalent to the emergent stellar luminosity) and 
the convective mach number (Mc)- Over the range of masses 
represented, Mc increases with luminosity (and mass), but 
only as ~ ^c^nv The yet weaker dependence of the con- 
vective turnover frequency on mass reflects the increase in 



Table 1. Core Convection Parameters 



^ Note that the simulations of iRoeers fc Glatzmaiej ||2006|) use 
an artificially high convective luminosity relative to the solar 
value. It is unclear whether this changes the functional form of 
the g-modc power vs. frequency; this needs to be studied in more 
detail. 

^ http://mesa.sourceforge.net/ 

■* We use a step function overshoot prescription, in which the 
convection zone is extended a distance of 30 per cent of a pressure 
scale height above the Ledoux boundary, with a constant diffusion 
coefficient. 



Mass'* Lconv' 

[Mq] [Lq] 



Mc 



[mHz] 



ZAMS TAMS 



2 


14 


4.37 X 10" 


'5 


0.05 


5.81 


1.58 


3 


75 


6.67 X 10" 


-5 


0.07 


5.91 


1.74 


4 


228 


8.75 X 10" 


-5 


0.08 


5.54 


1.71 


5 


579 


1.23 X 10" 


-4 


0.10 


4.73 


1.65 


6 


1.07 X 103 


1.40 X 10" 


4 


0.10 


4.63 


1.68 


7 


1.79 X 10^ 


1.57 X 10" 


4 


0.11 


4.57 


1.63 


8 


2.79 X 10^ 


1.73 X 10" 


4 


0.11 


4.12 


1.59 


9 


4.09 X 10^ 


1.88 X 10" 


4 


0.12 


4.22 


1.54 


10 


5.86 X 10^ 


2.07 X 10" 


4 


0.12 


4.19 


1.50 


12 


1.03 X lO" 


2.32 X 10" 


4 


0.13 


4.29 


1.43 


15 


1.98 X IQ-* 


2.59 X 10" 


4 


0.14 


4.43 


1.31 


20 


4.46 X 10-* 


3.00 X 10" 


4 


0.15 


4.35 


1.15 


25 


7.97 X 10-* 


3.28 X 10" 


-4 


0.15 


4.32 


0.99 


30 


1.23 X IQS 


3.46 X 10" 


4 


0.15 


4.36 


0.93 



" Initial stellar mass 

^ Convective luminosity on ZAMS 

^ Convective mach number on ZAMS 

Convective turnover frequency on ZAMS 
^ Minimum frequency for £ = 1 standing waves 

convective core radius with mass, since energy-bearing ed- 
dies span roughly the core radius. 

Fig. 13] shows propagation diagrams for our 2 and 10 
Mq ZAMS models (upper and lower panels, respectively). 
Each panel shows the Brunt- Vaisala frequency, A'^, in the 
radiatively stable envelope (blue, solid line), the Lamb fre- 
quency {Sj = £{£+ l)cl/r'^) for £ = 1,10 (green, solid 
line), and the outer convective turnover frequency in the 
core (red, dashed line). Note the ordinate units are linear 
frequency, i/. These diagrams highlight where g-modes and 
sound waves may propagate in the stellar interior, with g- 
modes propagatin g wherever a; < N,Se and sound waves 
where lj > N,Si (|Aerts et al.ll2010l : lUnno et al.lll989l ). The 
convective turnover frequency for all models is ~ 3 orders- 
of-magnitude less than the Brunt-Vaisala frequency in the 
radiative zone. Thus, modes with lu > luc propagate as 
high-order (short- wavelength) g-modes in the radiative zone. 
Lastly, the contrast between panels shows the growth of the 
convective core in fractional radius (and mass) with increas- 
ing stellar mass and the structural similarity among main 
sequence models with core convection zones. 



4 GRAVITY WAVE PHYSICS 

To identify potentially observable g-modes in massive stars, 
we must first determine whether waves excited at the 
convective-radiative boundary reach the stellar surface. 
Waves excited at the boundary travel through the radia- 
tive zone, where they are subject to damping by radiative 
diffusion. The radiative damping rate for a traveling g-mode 
is given by 

7rad(t^,^, r) =Krad(r) k^, 
16 a T{rf 



mesh_delta_coef f = 0.25 and varcontrol.target 



10" 



AiV 



(4) 
(5) 

(6) 
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Figure 1. Convective excitation spectrum vs. uj/uJc (abscissa) and I (ordinate), for g-modes in our 2 model. The spectrum shown in 
the left panel is for excitation by the convection zone, = 1, and the right for excitation dominated by an overshoot region having a width 
of 10 per cent of a pressure scale height, 77 = 0.1. Note the different scales on the abscissa in the two panels, since our parametrization of 
the power spectrum applies only above 10a;c for r) = 0.1 (see eon. [2ll. Also note values for log di?g/(d In oj din £) < are represented as 
white. The spectra for models with M > 10 Mq peak at higher £ values since r/H>2 at the edges of their convective cores; otherwise, 
spectra for different masses are qualitatively similar, with the overall energy injection rate scaled according to the value of Ale Lconv (see 
Tabled. 
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Figure 2. H-R Diagram showing the evolutionary tracks of our 
2 to 30 Mq model sequences spanning ZAMS to TAMS {Xc = 
0.01). The circles, triangles and squares mark the locations of 
models for which we have calculated convectively excited modes, 
and are in the evolutionary states given in the legend. Crosses 
mark every 1/lOth of the main sequence lifetime (i.e., the model 
age at TAMS). 



where A'^ = £{£ + I) and kr is the radial wavenumber. The 
last approximation (eqn. [6]) relies on the WKB dispersion 
relation for g-modes, which have lj <^ N, Si. 

We write a wa ve "optical depth" (as in 
iKumar fc Quataert|[l997l ) 



/r 



dr 



, 7rad(^,-g,r') 
^group 



where 

'group ~ ^ 



1^1 



AiV 



(7) 



(8) 




Figure 3. Propagation diagrams for 2 Mq (upper panel) and 
10 Mq (lower panel) ZAMS models. The solid blue line shows 
the Brunt- Vaisala frequency (N), the solid green lines the £ = 
1,10 Lamb frequency (Si), and the dashed red line shows the 
convective turnover frequency in the core {i^c}- 



Waves of a given frequency and degree deposit most 
of their energy at the radius where r^, ~ 1. As 7rad is 
inversely proportional to frequency, low frequency waves 
damp well inside the star while high frequency waves prop- 
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18 16 14 12 10 8 6 
logP [erg cm^^] 



Figure 4. Location, in log pressure units, where an g-mode of 
frequency u reaches unity optical depth due to radiative damping 
(see eqn.O or its outer turning point, for 2 Mq (upper panel) and 
10 M0 (lower panel) ZAMS models. Filled circles show the char- 
acteristic convective frequency in the convective core, diamonds 
mark frequencies with < 1 at the outer turning point (where 
they reflect and can set up standing waves), and the horizontal 
dashed lines give the minimum standing wave frequency for each 
Note that the outer turning point moves inward at high fre- 
quency as it is defined by where ui = Si. Colors in each panel 
correspond to as described in the legend in the lower panel. 



agate to the surface. The minimum frequency for standing 
g-modes satisfies the condition rw{ijJ,£, router) ^ 1, where 
router is the radius of the wave's outer turning point where 
uj = min (A'', Se). 

Fig. m shows the location where a wave of frequency 
f has an optical depth, Tu, ^ 1, for g-modes with i/ ^ Vc 
in our 2 and 10 Mq models. The horizontal dashed lines 
show the minimum frequency for standing g-modes, sepa- 
rating the low frequency waves that damp in the radiative 
envelope from those that reach their outer turning point 
and set up standing oscillations. Also shown are the char- 
acteristic turnover frequency in the convective core {vc, 
filled circles), and the outer turning points of g-modes with 
Tu,{u},£, router) < 1 (filled diamonds). At high frequencies, 
the outer turning point moves inward with increasing fre- 
quency since it is defined by where uj = St (see Fig. [3)l . 



5 QUASI-ADIABATIC STELLAR 
OSCILLATIONS 

For each model, we use the publicly available ADIPL^ adi- 
abatic stellar oscillation package to determine the eigen- 
functions for g-modes above the minimum frequency for 
standing waves. We then calculate the fiux perturbation and 
growth/damping rates, in the qu asi-adiabatic limit, a ccord- 
ing to the procedure described in IShiode et all l|2012l . eqns. 
1 - 10). 

Our quasi-adiabatic calculations are only valid in the 
regime where the thermal time of the overlying layers is 
greater than the mode period. Thus, we approximate surface 
perturbations by their value in the vicinity of the adiabatic 
cutoff, Rad, where the thermal time in the overlying layer is 
equal to the mode per iodlll Comparison with non-adiabatic 
growth/damping rates and eigenfunctions (Townsend & 
Teitler m prep) shows that our quasi-adiabatic calculations 
provide a good approximation for damping rates and surface 
amplitudes to within factors of a few. 

The short wavelength g-modes predominantly excited 
by convection are subject to large radiative damping. Those 
able to set up global standing waves are still linearly damped 
with lifetimes of ~ 1 — 10^ yr (with the exception of some lin- 
early driven modes). Damping rates are fit reasonably well 
by power laws in mode frequency, with indices of —3 to —8, 
where fiatter power laws apply for more evolved models. 
Deviations from the power-law arise primarily from opac- 
ity effects in the stellar envelope, whose relative importance 
depends on mode frequency. 

The convective excitation rate and mode lifetime to- 
gether set the equilibrium mode energy 

^-°-(-'^) = ^dl^N-(-'^)"^"' 

where En,„dc = a;V2 /y"'"^ [C? + A'Ch] dm, N,„(a;,^) is the 
number of modes in a logarithmic bin in u (at fixed I), 
and 7 is the damping rate of the mode as determined by 
all non-adiabatic effects (radiative diffusion, nuclear driving, 
and any convective viscosity). For g-modes with lj <^ N, 
Nn^{ij,£) ^ (n + //2)(2^ + 1), where n is the number of 
radial nodes in the mode eigenfunction l|Unno et al.lll989l ). 

In linear adiabatic stellar oscillation theory, the ampli- 
tude of the temperature and density perturbations, and the 
magnitude of the fluid displacement, are arbitrary. We con- 
vert the linear theory results into realistic predictions by 
normalizing the g-mode eigenfunctions using the estimate of 
the mode energy in equation [9] We then use these properly 
normalized eigenfunctions to calculate the disk integrated 
perturbations as follows. 



^ Distributed with MESA and separately at 
|http:/ /users-phys.au.dk/jcd/adipack.n/' 

' In detail, we take the median value of the perturbation in the 
range where this ratio is between 0.1 and 10. In our 10 Mq ZAMS 
model, the cutoff occurs at logP = 8.64, 8.46, 8.18, 7.73 ergcm"^ 
for our n = — 16, — 10, — 5, — 1 g-modes at frequencies of u = 
4.38,6.73,12.79,39.85 fiUz. 
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5.1 Surface perturbations 

For modes with 1^1, disk- averaging effects reduce the 
observable amplitude of oscillations as neighboring sur- 
face elements oscillate out of phase with one another. For 
each mode, we calculate the disk-integrated, limb-darkened, 
quasi-adiabatic perturbation to the bolometric magnitude 
and surface radial velocity for each stellar model (as in 
iDziembowskill 19771 ): 



h'jr + (26, - c) I 



\5m\ =1.087 

<5t;disk| = [Ul ^ + f£ Ch ^]r!^R ^ : 



(10) 
(11) 



where be,ce,ue, and vi are coefficients related to the effect of 
limb-darkening on the visibility of modes of a given £, 
are the radial and horizontal displacement eigenfunctions, 
and 5F/Fo is the fractional perturbation to the radiative flux 
due to the oscillation. We assume Eddington limb-darkening 
for simplicit y, and thus us e the b e,C£,ut, and ve coefficients 
tabulated in iDziembowskil l| 19771 ). 

We primarily use i5wdisk to compare our calculations 
with previous estimates and observational constraints on so- 
lar g-modes (see i]6.3p . For the more massive stellar models, 
which we show produce a much smaller velocity signal, we 
also consider a simpler estimate of the non-disk-averaged 
total velocity perturbation 



\Svtot \ = - [l^rUjf + l^hl^fY/^R^ 



(12) 



which might be connected to the non-thermal velocities 
needed to fit spectral line profiles in hot, massive stars: mi- 
croturbulence and macroturbulence. 

As a general rule, the procedure described above re- 
sults in a ratio of flux-to-velocity-perturbations of order 
Sm/Svdisk ~ 1 ^mag/ cm s~^ for £ = I modes. We find that 
this ratio varies between 0.1 — 10 ^mag/ cm s~^, with a trend 
towards larger values at lower frequencies. There is no clear 
trend with mass or stellar evolutionary state. 



6 RESULTS 

6.1 Brightness perturbations 

We find that, unless the excitation is dominated by an over- 
shoot region < 10 per cent of a pressure scale height thick, 
corresponding to a peak excitation frequency of > Wide, 
individual oscillation modes excited by convection produce 
brightness fiuctuations of Sm < 10 /imag (except at the 
TAMS, see gXT]). 

We thus focus on the effect of an ensemble of stochas- 
tically excited modes on the intrinsic variability of the star. 
This RMS magnitude perturbation is an incoherent sum over 
all modes in a given logarithmic bin in frequency: 



Srrir 



1 1/2 



J2 E [Sm{e,m,oj)Y 



A In (jj , all i m- 



(13) 



where the sum is over a logarithmic bin in uj and all I 
(though practically only ^ < 5 contribute). 

The left panels of Fig. [S] show the predicted 5m per- 
turbations due to individual 1=1 modes (top) and the 
total (5mrma (bottom) as a function of frequency for our 10 



Mq ZAMS models. The open circles in the upper left panel 
(for 5m) correspond to individual g-mode frequencies, while 
the circles in the bottom left panel are frequencies where 
5m^ras has been calculated from the individual modes. The 
flatter spectrum of the RMS magnitude perturbations, rel- 
ative to that of the individual t = 1 g-modes, results from 
summing the contributions from all modes in a given loga- 
rithmic frequency bin, including those of higher I (which are 
not plotted in the upper left panel; see equation I13|l . 

For ZAMS models, the combination of the convective 
excitation spectrum and that of radiative damping combine 
to produce a nearly flat spectrum of mode energies for the 
standing g-modes excited by convection, since E oc v~^^^^ 
and 7 oc . However, there is a strong decrease in surface 
amplitude with increasing frequency at fixed mode energy, 
which is well fit by 5m oc v~'^, due primarily to the outer 
turning point moving inward. This latter trend dominates 
the 5m spectra shown in the top panels of Fig.[Sl The spectra 
of flux perturbations is qualitatively similar to that shown 
in Fig. [5] for the range of stellar masses we investigated. 

The top panel of Fig. |S] shows the maximum (Jmmis in 
ZAMS models at each mass, with color corresponding to 
overshoot thickness and symbol size denoting the frequency 
of the peak flux perturbation. In more massive stars, the 
closer match between standing g-mode frequencies and the 
characteristic excitation frequency, and the overall larger lu- 
minosity and mach number, lead to larger brightness pertur- 
bations. However, the frequency of peak brightness, « 5 /iHz 
(~ 0.4 d^^), is roughly constant with increasing mass. 

6.1.1 Effect of stellar evolution 

Both the mach number in the core and the stellar luminosity 
increase during main sequence sequence evolution, leading 
to increased power input to g-modes and a larger charac- 
teristic convective frequency, Vc. In addition, the minimum 
frequency for standing waves decreases along the main se- 
quence as the thermal time at the mode outer turning point 
increases. Together these imply an increase in g-mode sur- 
face amplitudes with evolution, as well as the appearance of 
ever lower frequency modes. 

The top and bottom panels of Fig.[S]show the evolution 
along the main sequence of the power spectrum for individ- 
ual modes and iJmrms, respectively, in our 10 M© model. 
The RMS flux perturbation (lower panels) at fixed frequency 
increases with evolution, and the appearance of lower fre- 
quency modes at the surface leads to larger signals at lower 
frequencies than appeared on the ZAMS. 

The chemical composition gradient left behind by the 
receding convective core during the main sequence evolution 
produces an extended bump in the Brunt- Vaisala proflle; 
this results in enhanced trapping behavior near the end of 
the main sequence as modes separate into those with most of 
their energy in that composition-gradient region and those 
primarily in the radiative zone. This effect leads to the dif- 
ferences in the surface amplitudes of adjacent modes seen in 
the upper right panel of Fig. [S] 

As shown in Figs. [5] and [B] we predict g-mode power 
peaked at near 5/iHz on the ZAMS, but this peak ap- 
proaches ~ 1 /iHz at the end of the main sequence. However, 
in our evolved models, the minimum frequency for standing 
g-modes can be less than the minimum frequency for which 
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Figure 5. Individual £ = 1 g-mode (top panels) and RMS (lower panels; see eqn. I13I I surface brightness fluctuations as a function of 
frequency for convectively excited g- modes in our 10 Mq model at the ZAMS (left panels), midMS (middle panels), and TAMS (right 
panels). Individual 1 = 2 modes have similar amplitudes to those shown in the top panels for all phases, and £ = 3 modes approach 
parity at the TAMS. Open circles in the top panels mark individual mode frequencies, while in the bottom panels they show frequencies 
at which we have calculated (5mrnis. Note the decrease in the minimum standing wave frequency (i.e., the low frequency cutoff of the 
spectrum), and corresponding increase in the magnitude of the maximum brightness perturbation, with evolution. In our evolved models, 
the low frequency cutoff of the spectrum for r\ < 0.3 is set by the minimum excitation frequency for which our model is valid, Vcl'r), since 
this is larger than the minimum frequency for standing waves. There is likely to be convectively excited power at these frequencies (see, 
e.g., [Rogers &: Glatzmaier,,2006i) . but it is not accounted for here. 



our excitation model applies, uJc/ti, when the overshoot re- 
gion is < 30 per cent of a pressure scale height. This is 
shown clearly in the low-frequency cutoff of the spectra in 
Fig. O and in the lower panels of Fig. [6l where the mag- 
nitude of the maximum perturbation for 77 = 0.05,0.1 is 
smaller than that for larger r) and appears at a higher fre- 
quency. There is likely convective excitation power at these 
low frequ encies, but it is not accou nted for in our model 
(see, e.g.. [Rogers fc Glatzmaierll2006t ). 

If we assume the energy flux into modes below 
Vc/ri is constant, a s appe ars to be the case in the 
[Rogers fc Glatzmaiej (|2006l ) simulations, the equilibrium 
energy in modes below this cutoff scales as ~ i/*. However, a 
strong scaling of surface amplitude at fixed mode energy, as 
we find in the case of our 10 TAMS model, can lead to 
an increase in the maximum iJmrms by a factor of a few and 
a shift of the peak to slightly lower frequencies than appear 
in the lower panels of Figs. [S]and[Sl At yet lower frequen- 
cies, however, approaching the minimum for standing waves, 
fmin, the predicted flux perturbations decrease. 

For Kepler and CoRoT observations with micromag- 
nitude precision photometry, the RMS brightness fluctua- 
tions contributed by convectively excited g-modes should be 
detectable during the main sequence for all stars considered, 
if the driving is dominated by a convective overshoot region 
with width < 30 per cent of a pressure scale height at the top 
of the core. The expected amplitudes reach ~ 100s of /imag 
around z/ ~ 1 - 10 /^Hz (0.08 - 0.8 d"^). 



6.2 Velocity perturbations 

In addition to brightness perturbations, non-radial modes 
produce velocity perturbations which may be observed ei- 
ther in the disk-integrated radial velocity signature or via 
the motions they generat e where spectra l lines are formed 
( "microturbulence" ; e.g., ICantiello et al.l |2009| . and refer- 
ences therein). In massive main sequence stars, we are in- 
terested in the latter case; moreover, we are concerned with 
the aggregate effect of many modes, as discussed above for 
the brightness perturbations. 

In Fig. [T] we show the maximum predicted RMS ve- 
locity fluctuations as a function of initial stellar mass and 
evolutionary phase. For g-modes excited by the convective 
core, we predict surface velocity fluctuations that are always 
^ 1km s~^, whereas the observed non-thermal surface ve- 
locities in massive main sequence stars are > 1 — 10 km s~^. 
Even in our most optimistic scenario for excitation, with 
rj = 0.05, g-modes excited by core convection do not pro- 
duce sufficiently large velocities to explain the observed mi- 
croturbulence. 



6.3 Observability of solar g-modes 

We have also computed a solar model according to mod- 
ern asteroseismology constraints, in order to determine the 
surface velocities of solar g-modes implied by our excita- 
tion model. If the excitation is dominated by the highly- 
stratified solar convection zone, g-modes of a given frequency 
ij are primarily e xcited at the place in the convection zone 
where a; ~ (r) (|Kumar et al.lll999l ). Thus high frequency 
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Figure 6. Peak RMS surface brightness fluctuations (5mrms for 
convectively excited g-modes, calculated according to eqn. 1131 as 
a function of mass. Color corresponds to the value of r) and size 
to mode frequency, as described in legend. Panels show results 
for models at their labeled evolutionary states. In the evolved 
models, the minimum frequency for excitation at small r) {vc/r)) 
is larger than the minimum frequency for standing waves (t'niin; 
Table 1), leading to a suppression in the total flux perturbation 
for the smallest values of r). 
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Figure 8. Disk-integrated velocity perturbations for £ = 1 solar 
g-modes. For comparison, £ = 2, 3 modes have smaller velocity 
perturbations by factors of ~ 3 and 10, respectively. The blue 
points show the amplitudes of modes assuming that the highly- 
stratified solar convection zone dominates the exciation, as com- 
pared to the green and red points which give the result of excita- 
tion by an overshoot region of 30 or 5 per cent of a pressure scale 
height, respectively. 



g-modes are predominantly excited near the solar surface, 
while low frequency g-modes are excited near the radiative- 
convective boundary in the interior. Assuming the convec- 
tion zone behaves like a polytrope with index n « 1.5, equa- 
tion [1] still applies, but with ^max = r / H (lu/u]c)'^^^ , where 
r and H are evaluated at the interior radiative-convective 
boundary where g-modes begin to propagate, a = 7/2, 6 = 2, 
and rj — 1. 

However, if the excitation is instead predominantly due 
to a thin overshoot region of width rjH at the base of 
the convection zone, then we assume the excitation spec- 
trum at fixed position {a = 13/2), with < 1. In Fig. 
[S] we show the disk-integrated g-mode velocity amphtudes 
that result from using the convection-zone-dominated and 
overshoot-dominated excitation spectra. In the overshoot 
case, we choose values rj = 0.05, 0.3 that span the range sug- 
gested by helioseismic investig ations and numerical simula- 
tions of solar convection (e.g., Christensen-Dalsgaard et al.l 



I2OIII : [Rogers fc Glatzmaieij|2005l . boosf ) 



In all reasonable excitation scenarios, we find ampli- 
tudes of 5«disk 5, 0.3mm s~^, which are consistent with 



those derived by Kumar et al.l (|l996l ). Our results are also 
consistent with the observational upp er limits that place 
^fdisk ^ 10 mm s"'^ at ^ 100 /iHz (|Appourchaux et al.l 
[20l3). 

The intensity-to-velocity ratio we find for these high 
frequency solar g-modes is Sm/5v < 1 /xmag/ cm s~^, so 
that typical intensity perturbations are < 10~^ /^mag. 
This is also well below the published upper limits from 
lAppourchaux et al] l|2000l ) of ^ 0.5 ^mag. 



Figure 7. Peak RMS surface velocity fluctuations convectively 
excited g-modes as a function of mass. Symbols and colors are as 
in Fig.E] 



7 DISCUSSION 

Using a simple parametrization for the convective excita- 
tion of gravity modes, we predict that the convective cores 
of stars with masses > 2 Mq will excite observable surface 
brightness fiuctuations while on the main sequence. If the ex- 
citation is dominated by an overshoot region having a width 
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of < 30 per cent of the pressure scale height at the top of the 
convective core, these stars wiU exhibit intrinsic photometric 
variability with amplitudes up to 10s of micromagnitudes at 
frequencies of 5 - 10 /xHz (0.4 - 0.8 d"^) on, or just red of 
the solar metallicity ZAMS. These flux variations are largest 
for massive stars with M > 5 Mq (Fig. [B]). 

As the stars evolve along the main sequence, we predict 
a strong increase in the flux variability, reaching 100s of 
micromagnitudes at frequencies < 10 fiRz, along with the 
appearance of power at ever lower frequencies, down to ~ 
IfiUz (0.08 d"^) at the TAMS. We predict no observable 
accompanying velocity signal. 

Given the relatively low amplitudes of individual g- 
modes, we have framed our predictions in terms of the RMS 
flux variations at a given frequency, rather than the am- 
plitudes of individual stellar normal modes. The latter are 
smaller by factors of ~ 3 — 30 at these frequencies (with 
larger corrections at lower frequencies), making it challeng- 
ing to detect individual normal modes. However, as the du- 
ration of observation for Kepler and CoRoT stars increases, 
it may become possible to detect individual normal modes 
in a periodogram analysis, since the long lifetimes of these 
modes ensures that they will be come more promin ent in 
longer duration observations fe.g. iDupret et allbOOgl ). 

Our excitation model agrees with previous predic- 
tions that the surface amplitudes of solar g-modes are < 
0.3 mm s^^ (Fig. |8]), thus implying that photometry of mas- 
sive mam sequence stars provides the best potential window 
into the convective excitation of g-modes. 



7.1 Theoretical uncertainties 

7.1.1 Overshoot excitation 

Detailed asteroseismic modeling provides a best fit width for 
the solar overshoot region ( "tachocline" ) of between 5 and 
40 per cent of a pressure scale height at the base of the con- 
vection zone, depending on the ove rshoot prescription used 
l|Christensen-Dalsgaard et al.l 20 111, and references therein) . 
Numerical simulations bv lRogers fc Glatzmaieil l|2005l . l200a ) 
imply an overshoot thickness at the small end of this range 
^ 0.05//. However, asteroseismic results from other solar- 
like stars suggest that the thickness of the overshoot layer 
may not be universal eve n among similar stellar types (e.g., 
iLebreton fc Goupilll2012l '). 

For main sequence stars with core convection zones, nu - 
merical simulations by [Browning. Brun. fc Toomrd (|2004l ) 
find that the overshoot region has, in our terms, r/ < 0.2. 
Results from asteroseismic modeling of massive main se- 
quence pulsators are broadly consistent with overshoot of a 
few tenths of a pressure scale height, but show a r ange of best 
fit pa rameters from consistent with zero (e.g., lAerts et al.l 
I2OIII ) to more than 40 p er cent of a pressure scale height 
(e.g.. iBriguet et aTlbOOTl ): furthe r corroboration is found in 
isochrone fitting of open clusters l|Brott et al.ll201ll . and ref- 
erences therein). This broad base of observational results, 
in combination with arguments presented here, suggest that 
g-mode excitation by convective overshoot could produce 
observable surface brightness perturbations in massive main 
sequence stars. 



7.1.2 Effects of rotation on g-modes 

Usin g typical observed surface velocities for massive stars 
(e.g. I Wolff et al.]|2006l ). the corresponding rotation frequen- 
cies are given by 



SI ~ 3.5 ^Hz 



150 km s-i 



R 

R0 



(14) 



i.e., S7 ~ 3 (1) fiUz for M = 2 (10)Mq. 

The effect of the Coriolis force on g-mode propagation 
can be seen in th e WKB dispersion relation (as in e.g., 
iKumar et ahll 19991 . eqn. 21): 



2 ,(k-n)^ 

cj^ ^ 4 ^ ' + N 



fe2 



(15) 



Thus, g-modes only propagate where cj > 2(n • fc), and 
modes with frequencies a; < are confined to the equator. 

As discussed in ^ much of the convectively excited 
power is in modes with frequencies of 1 — 5 ^Hz, com- 
parable to the median rotation frequencies of massive main 
sequence stars. Thus we expect that rotation will have a 
non-negligible effect on the surface manifestation of convec- 
tively excited g-modes (it may also change the excitation 
of g-modes since rapidly rotating convection has different 
statistical properties than non-rotating convection). This is 
particularly true as the star evolves to the TAMS. Given the 
potential importance of rotation, it would be particularly in- 
teresting to compare high precision photometry of slowly vs. 
rapidly rotating main sequence A-0 stars. 

Additionally, if the surface of the star is rotating more 
rapidly than the core, the star may possess critical layers 
that significantly increase the radiative damping of outgoing 
g-modes and preclude the establishment of standing waves. 
Such critical layers may in fact be set up by the angu- 
lar momentum carried by the g-modes themselves (see e.g.. 



lar momentum carried by 
[Rogers, Lin, fc Laullioia ) 



7.1.3 Effects of stellar model parameters 

To assess the effects of stellar modeling parameters on our 
predictions, we evolved a 10 Mq star with variable mix- 
ing parameters and rotation. We tested mixing length pa- 
rameters aMLT ~ 1,1.5,2; boundaries determined by the 
Schwarzschild and Ledoux criteria without overshoot; and 
the Schwarzschild criterion with 30 per cent overshoot and 
initial surface rotation velocities of 50, 100, 200km s~^ (with 
and without the chemical mixing induced by Taylor- Spruit 
magnetic fields). 

We find that the uncertainty in our predicted mode 
amplitudes at fixed central hydrogen mass fraction due to 
stellar evolutionary parameters may be up to a factor of 
~ 5, with the models presented in our figures roughly cor- 
responding to the median of the distribution. Changing 
Qmlt produces the largest effect, as larger auur corre- 
sponds to smaller ujc and larger A4c', the overall normal- 
ization of the convective excitation spectrum scales as Q^lt 
for a = 13/2 (see eqn.[T}. Convective boundary definitions 
can also change the ratio of Wmin/wc by up to a few percent 
on the ZAMS, changing the fraction of the convectively ex- 
cited g-mode power that reaches the stellar surface by < 25 
per cent. By the end of the main-sequence these effects can 
amount to an order of magnitude dispersion among models. 
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as the determination of convective boundaries affects the 
progression of stellar evolution. 



7.1.4 Turbulent power spectrum 

In their study of th e detectability of solar g-modes, 
iBelkacem et al.l l|2009al ) employ a Lorentzian eddy-time- 
correlation function for coupling the turbulent convective 
motions to g-mode excitation. They find that this results in 
a 30-fold increase in solar g-mode amplitudes, as compared 
to t he assumptions used in our analysis (and previous work; 
e.g.. lKumar et al]l 19961 ). ISamadi et al.l (|2010l ) also employed 
this formalism to study convectively excited g-modes in 10, 
15, and 20 Mq main sequence models having Xc = 0.5. They 
find amplitudes of ~ 10 /imag for individual 1=1 g-modes 
in these stars, also ~ 30 times larger than the amplitudes in 
our 77 = 1 case (see our Fig. [5] and their fig. 2). 



7.1.5 Near-surface convection zones 

Massive main sequence stars of near solar metallicity have 
vigorous ne ar-surface convection zones due to the iron opac- 
ity bump jCantiello et al.l l2009h . These convection zones 
are another potentially important source of g-modes near 
the stellar surface, having mach numbers of > 0.01. Near- 
surface convection zones have characteristic frequencies of 
10s to 100s of /iHz, but the small scale-heights near the 
stellar surface imply that nearly all of the convective power 
is at larger characteristic £ > 30. These high £ g-modes are 
unlikely to produce any signific ant surface brightness per- 
turbations, but, as suggested bv lCantiello et al](|2009l ). the 
velocity field associated with these waves might account for 
the micro/macro-turbulent velocity fields inferred via spec- 
troscopic modeling. We leave a detailed investigation of the 
excitation and propagation of these modes to future work. 



7.2 Observational prospects 

7.2.1 Distinguishing stochastically excited modes 

Some of the stars we have investigated lie in the SPB in- 
stability strip, where stars are observed to pulsate in low- 
degree g-modes with frequencies of 0.6 — 3d~^. These lin- 
early excited modes have frequencies very near the peak of 
the stochastically excited modes. It is thus important to de- 
termine how to observationally distinguish between linearly 
and stochastically excited modes. 

Unlike p-modes stochastically excited by envelope con- 
vection zones (i.e., solar- like oscillations), the convectively 
excited modes explored in this work have linear damping 
timescales (primarily due to radiative diffusion) of years or 
longer, much longer than a typical observation timescale. 
Thus the spectral broadening and time-frequency diagram 
characteristics that distinguish stochastically-driv en modes 
in the former case cannot b e leveraged here (e.g.. ISeddinel 
l201ll : lBelkacem et al.ll2009bl '). We must instead rely on the 
fact that convectively excited g-modes should be observable 
in stars outside the linear instability strips and at frequen- 
cies that should otherwise be damped in linear analyses. 



7.2.2 Have convectively excited g-modes already been 
detected? 



lAerts et all (|2009l ') have shown that the "macroturbulent" 
velocity fields observed in spectra of blue supergiants can be 
explained by the collective effect of an ensemble of excited 
g-modes, which span a similar frequency range to that pre- 
sented here. More recent ongoing work bv lSimon-Dfaz et al.l 
I2OIII ), has shown that macroturbulence is present in O and 
B dwarf stars as well, suggesting that perhaps the same 
low-frequency pulsations are present on the main sequence. 
We have shown however, that if individual mode ampli- 
tudes are set by the competition between convective exci- 
tation by the core and radiative damping in the envelope 
(see eqn.[9|, the aggregate effect is too small to explain the 
observed macroturbulence (see Fig. [T]). The source of these 
turbulent motions could instead be g-modes excited by near- 
surface convect i on zo nes, as discussed in i|7.1.5l above and 
ICantiello et al.l (|2009l ) . Another interesting possible source 
is the photon bubble in stabilities that occu r in the atmo- 
spheres of massive stars l|Turner et al.|[20o3 ). 

There is tantalizing evidence for convectively excited 
g-modes in recent observation s of A, B and O stars 
with Kepler (iBalona et al.ll201ll : lUvtterhoeven et al.ll201ll : 

ul ' 



Blomme et al.l 



20lj) 



In their study of B-star pulsators. 



Balona et alT(|201ll ) find evidence for theoretically unex- 



pected, low frequency pulsations in SPB stars and SPB//3 
Cep hybrid pulsators. Some of these stars lie outside the 
classical instability strips for low-order modes, potentially 
implying that a previously unexplored driving mechanism is 
at work in these stars. 

In studies of O-stars, iBlomme et al.l l|201ll ) have found 
an unexpected "red-noise" component to the stellar photo- 
metric power spectrum, even after instrumental corrections 
have been applied. This noise appears at the frequencies and 
amplitudes that agree with our predictions; however, the in- 
ferred mode lifetimes of hours to days for the "red-noise" 
modes are incommensurate with those of the g-modes ex- 
plored here, which have lifetimes of years to Myrs. 

As a general comment, we note that many studies re- 
move or disregard power at low-frequencies below 0.5 or 
0.2 d~^, depending on the study. They appeal to instrumen- 
tal effects both known and un-characterised. However, we 
urge caution: the complete removal of this low frequency 
power may take with it the signatures of convectively ex- 
cited g-modes. 
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